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Abstract. We consider inflaton couplings to MSSM flat directions and the 
thermalization of the inflaton decay products, taking into account gauge symmetry 
breaking due to flat direction condensates. We then search for a suitable curvaton 
candidate among the flat directions, requiring an early thermally induced start for 
the flat direction oscillations to facilitate the necessary curvaton energy density 
dominance. We demonstrate that the supersymmetry breaking A-term is crucial 
for achieving a successful curvaton scenario. Among the many possible candidates, 
we identify the uidd flat direction as a viable MSSM curvaton. 
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1. Introduction 

Minimal Supersymmetric Standard Model (MSSM) provides nearly 300 gauge in- 
variant flat directions, whose potentials are vanishing in a perfect supersymmetric 
(SUSY) limit (for a review, see PP). However, in the early Universe SUSY is broken 
and the flatness of the potentials is spoiled, but the directions still remain flat as 
compared to the curvature scalar of the Universe. Soft SUSY breaking terms induce 
mo ^ ©(TeV) masses to the flat direction. Flatness is also lifted by non-renormalizable 
superpotential terms, the form of which is dictated by the gauge properties of the flat 
direction. When the Hubble expansion rate becomes equal to the low energy SUSY 
scale, all the flat directions are trapped in their local minima, after which they start 
to oscillate and ultimately decay. The flat directions may a play a crucial role in 
cosmological issues ranging from baryogenesis to dark generation matter. Recently 
it has also been pointed out that they are also important for understanding the full 
thermalization history of the Universe after inflation |21 EI ■ 

An MSSM flat direction may also account for generating adiabatic density 
perturbations ^ [H |7] through the curvaton mechanism |S1 |H1 EI ■ Since during 
inflation the Hubble expansion rate Hj ^ toq it does not cost anything in energy, 
quantum fluctuations are free to accumulate (in a coherent state) along a flat direction 
and form a condensate with a large VEV, (pQ. Because inflation smoothes out all the 
gradients, only the homogeneous condensate mode survives. However, the zero point 
fluctuations of the condensate impart a small, and in inflationary models a calculable, 
spectrum of perturbations on the condensate. 

After inflation, H (x t^^ and the flat direction stays at a relatively large 
VEV due to large Hubble friction term imtil H ~ mo, whence the flat direction 
condensate starts oscillating around the origin with an initial amplitude ipQ. From 
then on \ip\ is redshifted by the Hubble expansion oc H for matter dominated and 
^ ^3/4 j-adiation dominated Universe. The energy of the oscillating flat direction 
may eventually start to dominate over the inflaton decay products. When the 
flat direction decays, its isocurvature perturbations will be converted to the usual 
adiabatic perturbations of the decay products, which thus should ultimately contain 
also Standard Model (SM) degrees of freedom. However, such an evolution is not 
automatic but is subject to various constraints 01 Ej, with the outcome that in general 
it is very difhcult to have a MSSM flat direction curvaton that could give rise to the 
desired spectrum of density perturbations. 

The key ingredient in these considerations is the process of thermalization of the 
ambient plasma. The state of the plasma before complete thermalization depends on 
the nature of flat direction which is developing a large VEV. For instance, MSSM 
Higgses developing a VEV give masses to the gauge bosons of S'C/(2)w x U{1)y, but 
gluons remain massless. The decay products can then reach partial thermalization 
through the gluons. Therefore one needs to study various possibilities case by case. 
It is important to know the details of the inflaton couplings to the MSSM degrees as 
well, since not all the couplings are renormalizable. 

Thermal effects are important for the MSSM curvaton mechanism I11[I12| because 
large thermal corrections can evaporate or dissociate the MSSM curvaton candidate 
before it has a chance to dominate the energy density. These effects obviously depend 
on the thermalization rate. Recently it has been pointed out that the VEVs of 
the gauge invariant flat directions actually slow down the thermalization rate of the 
inflaton decay products |51|3I. This is due to the fact that, in general, the flat direction 
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VEV gives masses to gauge bosons and gauginos which decrease the rates governed 
by 2 <-^- 2 and 2 — > 3 scatterings mediated via gauge bosons and gauginos. After the 
infiaton decay | the initial plasma will be in a quasi-thermal phase for a long period 
during which kinetic equilibrium is established but not the full chemical equilibrium 
and the initial plasma is far from thermal equilibrium. 

Another consequence of thermal corrections is that if the flat direction survives 
long enough, then thermally induced mass corrections can trigger early oscillations, 
which may hasten the dominance of the curvaton over the infiaton decayed products. 
These subtleties has never been systematically considered in the context of an MSSM 
curvaton. It is the purpose of the present paper to make use of the results of ref. 12 13] 
and study these effects properly. 

In previous work the evolution of the MSSM curvaton has been discussed 0] E] 
without due consideration of the SUSY breaking A-term in the potential, which is 
induced by the presence of the non-renormalizable terms in the superpotential. Here 
we take the A-terms into account and show that they greatly modify the dynamical 
evolution and in particular the initial amplitude of the curvaton oscillations. As a 
consequence, many of the obstacles for an MSSM flat direction curvaton are now 
removed, and considering the various suitable directions case by case we are able to 
pick out the Uidd flat direction as the most promising curvaton candidate. 

The paper is organized as follows. In Section|21we discuss how the infiaton couples 
to the MSSM degrees of freedom, emphasizing that there are certain combinations of 
slepton and squark superfields that do not have renormalizable couplings to a gauge 
singlet infiaton. In Section |31 we disseminate the potential flat direction curvaton 
candidates and choose those which give rise to thermal effects that can trigger early 
curvaton oscillations. Section ^ explains the importance of the A-teim for MSSM 
curvaton dynamics, while in Section ]E\ we discuss the thermal corrections to flat 
directions, which singles out the udd direction as the curvaton. In Section|Hlwe discuss 
the dynamics of the udd curvaton and find its initial amplitude and a lower bound 
on the renormalizable infiaton coupling. Our conclusions are presented in Section 

2. Infiaton decay 

In almost all F-term or I?-term models of inflation the inflaton, $, is considered to be 
an absolute gauge singlet §. Then the pertinent question is how the infiaton couples 
to the matter. 

Gauge symmetry implies that the infiaton must be coupled to a gauge-invariant 
combination of the SM fields and their SUSY partners. The field content of MSSM is 
governed by the following superpotential: 

Wmssm = A„QH„u + AdQHdd + AeLHde + /iH„Hd , (1) 

where Hu, Hd, Q, L, u, d, e in Eq. are chiral superfields representing the two Higgs 
fields (and their Higgsino partners), LH (s)quark doublets, RH up- and down- type 
(s)quarks, LH (s)lepton doublets and RH (s)leptons respectively. The dimensionless 

J Within SUSY the inflaton predominantly decays perturbatively. Non-perturbative decay of the 
inflaton to MS(SM) or through selfcoupling 1131 via parametric resonance is kinematically blocked 
because the inflaton decay products can couple (in a gauge invariant fashion) to MSSM flat directions. 
Large VEV of the flat direction induces large masses to the inflaton decay products which blocks 
preheating completely 1141 . 

§ One exceptional example is ref. 1151 . where the gauge invariant combinations of SO{N) flat directions 
are responsible for driving inflation. 
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Yukawa couplings A^, A^j, Ae are 3x3 matrices in the flavor space, and we have omitted 
the gauge and flavor indices. The last term is the ^ term, which is a supersymmetric 
version of the SM Higgs boson mass. 

There exist two gauge-invariant combinations of only two superfields: 

HuHd , HuL. (2) 

The combinations which include three superfields are: 

HuQu, HdQd, HdLe , QLd, udd, LLe. (3) 

SUSY together with gauge symmetry requires that the inflaton superfield $ is coupled 
to these combinations or to their superpositions. The terms $HuH<j and $HuL have 
dimension four, and hence are renormalizable. On the other hand, the interaction 
terms that couple the inflaton to the combinations in Eq. Q have dimension five and 
are non-renormalizable. They can arise after integrating out the heavy degrees of 
freedom which are associated to the new physics at high scales (for example GUT). 
Such terms are also generically induced by gravity, notably in supergravity, in which 
case they will be Planck mass suppressed 1 1 . 

The simplest case arises when the inflaton is coupled to matter via superpotential 
terms of the form: 

;i*HuHd,/i*Hu*, (4) 



where we have defined 

* = i ^ hiU , h 



1/2 



(5) 



Here 1 < i < 3 is a flavor index and h can be as large as 0(1). 

It is important to realize that one has to go from the {Li, L2, L3} basis, in which 
Ae in Eq. are diagonalized, to the {'4', L'l, L'2}, where L'l and L'2 are orthogonal 
to '3/, and hence have no renormalizable couplings to the inflaton; only ^ couples to 
the inflaton. (for a discussion on multi-flat directions, see 16 ). This is just a unitary 
transformation. 

Besides the SM gauge group the MSSM Lagrangian is also invariant under R- 

O D _|_ T _|_0 O 

parity (—1) . Preserving R-parity (at least) at the renormalizable level further 

constrains inflaton couplings to matter. Note that H„Hd is assigned +1 under R- 
parity, while has the opposite assignment —1. Therefore only one of the couplings 
in Eq. preserves R-parity: €>HuHrf if the inflaton superfield carries no lepton or 
baryon number, which is generically the case. In contrast, the 4>Hu'4' term could be 
relevant for e.g. models where the RH sneutrino plays the role of the inflaton |17j . 

In the context of a curvaton scenario the fact that the inflaton can have 
renormalizable couplings to (some of) the MSSM fields is an auspicious sign. The 
interactions in Eq. result in an inflaton decay rate: 

On the other hand, non-renormalizable interactions from Eq. Q result in a much 
smaller decay rate ~ m^/Afp. Such a late decay would imply a longer period of 

II It is possible that the inflaton decays mainly to another singlet (for example, the right-handed 
neutrino) suporflold. Then, since the ultimate goal is to create matter, this singlet must be coupled 
to MSSM flelds. Supersymmetry and gauge symmetry again require couplings to gauge-invariant 
combinations in Eqs. and along the lines above. 
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inflaton domination during which the ratio of the curvaton energy density to the total 
energy density of the Universe does not grow. Moreover, a small decay rate would 
result in a low reheat temperature, and hence would diminish the thermal effects 
necessary to trigger early oscillations of the curvaton field, as will be discussed below. 

3. Curvaton candidates 

Having discussed the inflaton couplings to matter, let us now attempt to identify the 
curvaton candidates within MSSM. An important point is that generating density 
perturbations of the correct size from curvaton decay requires that during inflation 
its VEV (fo ^ lO^iJi. Here Hi is the Hubble expansion rate during the inflationary 
epoch. 

A curvaton is a light modulus, lighter than the Hubble expansion rate during 
inflation. The total potential during inflation is given by 

Vtotal = V{I) + Vi^) (7) 

where V{I) is due to inflation and V{(p) is due to the curvaton. The curvaton carries 
isocurvature perturbations which sources the curvature perturbations. In order not 
to have any residual isocurvature perturbations left over the curvaton must decay 
dominantly into the SM degrees of freedom. The interesting quantity to study is the 
ratio of the perturbation and the background field value of the curvaton, since this 
is related to the curvature perturbation [S] El , which can be constrained from the 
amplitude of the CMB anisotropy. If the perturbations in the curvaton do not damp 
then the final curvature perturbation will be given by 

S^^^ 10-5 , (8) 

finf 

where 10"^ arises from the COBE normalization. In this we will present a scenario 
where the perturbations are not damped away and the curvaton dominates while 
decaying. 

There are three important points to note: 

• The flat direction that plays the role of the curvaton cannot include those fields 
which have a renormalizable coupling to the inflaton. The reason is that such fields 
generically acquire a mass h(j) 3> -ffm/ through their couplings to the inflaton, and 
hence cannot develop the required VEV. Therefore Eq. ^ implies that a curvaton 
candidate cannot include H^, Hd and ^I'. As pointed out earlier, such a late decay 
would be harmful to the curvaton scenario. 

• A curvaton candidate should not induce a mass > for the inflaton decay 
products, otherwise the two-body inflaton decay will be kinematically blocked. The 
decay will be delayed until the relevant flat direction has started its oscillation and 
its VEV has been redshifted to sufficiently small values 

• The flat direction which plays the role of the curvaton should not break all of the 
SM gauge symmetry. The reason is that in this case all gauge integrations will 
get decoupled and full (i.e. both kinetic and chemical) thermal equilibrium will be 
established much later than the inflaton decay '5/ %. However thermal effects 
which are necessary to trigger early oscillations of the flat direction \1^ 

^ If only part of the SM group is broken, gauge interactions of the unbroken part will bring fields 
which carry those gauge quantum numbers, and subsequently other fields, into full equilibrium. 
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crucially depend on the presence of a thermal bath consisting of some MSSM 
degrees of freedom in full equilibrium when the infiaton decay is completed. 

Note that a flat direction, which is a combination of the squark and slepton fields, 
has couphngs to other MSSM fields through which it induces a large mass for them. 
D-terms lead to flat direction couplings of gauge strength to the gauge flelds and 
gauginos associated with the gauge (sub)group which is spontaneously broken by the 
flat direction VEV, as well as to the orthogonal directions and their supersymmetric 
partners. F-terms result in flat direction couplings of Yukawa strength to the scalars 
that are not included in the monomial representing the flat direction, as well as to 
their fermionic partners. 

To elucidate these points, let us consider the particularly simple example of the 
H^Hd, flat direction+. In this case the flat direction and the orthogonal directions are 
defined as + ffd,2) /\/2, and {Hu^i - Hd,2) / V^, Hu,2 and H^.i, respectively, 

where 1,2 are the two components of iJ„ and iJ^. The flat direction breaks 
SU{2)w X U{1)y down to [/(l)cm- Then the gauge fields and gauginos associated 
with the broken subgroup, as well as (i?ti,i — Hd,2) /V^, Hu,2, Hd.i and their fermionic 
partner, acquire a mass ~ gi^p), where g is a gauge coupling. The Qi multiplets, which 

are coupled to both i7„ and Hd, acquire a mass ^(A^ ^ + ^){ip) while u^, and 

multiplets obtain respectively masses A„^i((p), Xd,i{(p), and Ae,,;(iy9). 
Let us denote the flat direction superfield by ip. Then we have 



where '^1,2,3 are MSSM superfields. In general the relationship < Hi„f holds *. 
For the curvaton mechanism to work, we require (fmf ^ lO^Hinf- Note that the VEV 
of the fiat direction induces VEV dependent SUSY preserving masses to the MSSM 
particles. Therefore, for the inflaton decay to be kinematically allowed, one needs 



if the #H„^' coupling is allowed. 

The above conditions considerably restrict the curvaton candidates within the 
MSSM as only the flrst generations of (s)quarks and (s)leptons which have a Yukawa 
coupling <^ 10"^. 

Let us now identify the fiat directions which satisfy all of the above mentioned 
requirements. An important point is that the acceptable fiat directions should include 
only one Q or one u. The reason is that D- and F-flatness of directions which involve 
two or more Q and/or u requires them to be of different flavors (for details, see [18\). 
This implies the presence of up- type squarks from the second and/or third generation 
which, according to Eq. @, will lead to Ai > 10"'^. This violates the condition for two- 
body infiaton decay via either of the $H„Hd or #H„^' terms, given in Eqs. (|10lll|) . A 
large number of MSSM fiat directions will therefore be excluded by this consideration. 
The only fiat directions with Ai < 10~^ are as follows: 

+ Note that, as we mentioned, this flat direction cannot obtain a large VEV. We only consider this 
example to demonstrate how the flat direction VEV induces mass to other MSSM flelds. 
* This is strictly true in models of chaotic inflation and (supersymmetric) hybrid inflation. It is 
possible to have ^ ^^in/ ™ models of new inflation, but this is a rather contrived situation. 



W^MSSM D AiH„<^Si + A2Hd<pE2 -I- Aa^yjSa, 



(9) 



Ai,A2 < 10-^ 
if the #H„H(i coupling is allowed by R-parity, and 
Ai,A3 < 10-^ 



(10) 



(11) 
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• udd : This monomial represents a subspace of complex dimension 6 I18j . D- 
flatness requires that the two d are from different generations (hence at least one 
of them will be from the second or third generation). This implies that A2 > 10"'^, 
see Eq. (O, for all flat directions classified by udd. As a consequence the two- 
body inflaton decay via the ^HuJld term will be kinematically forbidden. Note 
however that these flat directions are not coupled to Li_2.3- Therefore the inflaton 
decay via the €>Hu4' term can proceed for the Uidd directions, with ui being the 
RH up squark. We also note that these directions leave the S'J7(2)w unbroken, 
so that the S'?7(2)w degrees of freedom can completely thermalize. 

• QLd: This monomial represents a subspace of complex dimension 19 US]. F- 
flatness requires that Q and d belong to different generations. Then, since Q and 
d are both coupled to H^, Eq. © implies that A2 > 10~^. Thus two-body inflaton 
decay via the #H„Hc( term will be kinematically forbidden. On the other hand 
it can proceed via the $H„* term for QiL'd directions, where Qi is the doublet 
containing LH up and down squarks. Note that here we have to rotate to the 
{'S', 'L'1,1'2} basis where only L'^ and L'2 can acquire a large VEV during inflation. 
We also note that these directions completely break the SU{2)^^ x J7(1)y, but 
leave a SU{2) subgroup of the S'[/(3)c unbroken. Therefore the associated color 
degrees of freedom can completely thermalize. 

• LLe This monomial represents a subspace of complex dimension three |18j . D- 
flatness requires that the two Ls are from different generations, while F-flatness 
requires that e belongs to the third generations (therefore all the three lepton 
generations will be involved). Eq. (|5J) then implies that A2 — 10~^, which 
kinematically blocks two-body inflaton decay via the $H„H(; term. The decay 
can nevertheless proceed via the $H„^ term. However, since 5* cannot develop 
a large VEV, we should actually rotate to the {^',L'i,L'2} basis instead and 
consider the L'iL'2e monomial. Such a change of basis has no impact on D- 
flatness but will affect the F-flatness. The reason is that in general Ae is not 
diagonal in the {'4', L'l, L'2} basis, and hence L'^ and will be coupled to all 
flavors of e. This can be circumvented if the inflaton dominantly couples to one 
of the Li and ^' is mainly L^. More specifically, a feasible curvaton candidate 
will be obtained along the L2L3ei direction if SI/ Li. For this flat direction we 
have Ai = and A3 ~ 10^''' (see Eq. ©). This implies that two-body inflaton 
decay will proceed via the $H„Li term without trouble. We note that this flat 
direction completely breaks the electroweak symmetry SU{2)^ x [/(1)y, while 
not affecting 5[/(3)c- Therefore color degrees of freedom can fully thermalize. 
In passing we note that this monomial could be useful to generate primordial 
magnetic field, see |19j . 

• LLddd : This monomial represents a subspace of complex dimension three [THl . 
D-fiatness requires that the two Ls and the three ds are all from different 
generations. This implies that that A2 — 10^^ and Ai = A3 = (see Eq. Q). 
Therefore two-body inflaton decay via the ^HuH-d term will be kinematically 
forbidden, but the decay can proceed via the $H„'S' term. Again note that we 
should rotate to the {^,L'i,L'2} basis, and hence the L'iL'2ddd direction will 
be the relevant direction. However this direction breaks all of the SM gauge 
group. This results in late thcrmalization of the Universe [H] and the absence 
of thermal effects which, as a consequence, does not yield early oscillations of 
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the flat direction. This is an undesirable feature for our scenario in which, the 
curvaton oscillations are triggered by thermal effects, as we will discuss in the 
next Section. 

To summarize, after taking all considerations into account, we find that the most 
promising candidates for the curvaton within MSSM are the Uidd and QiL'd (and 
possibly L2L3ei) flat directions. 

4. Importance of the A-term 

Let us now discuss the role of the ^-term in constructing a curvaton model. In models 
with gravity and anomaly mediation, SUSY breaking results in the usual soft term, 
mo|(/3|^, in the scalar potential with mo ~ 100 GeV — 1 TeV. There is also a new 
contribution arising from integrating out heavy modes beyond the scale M, which 
usually induces non-renormalizable superpotential terms of the form 

where tp denotes the superfleld which comprises the flat direction ip. In general M could 
be the string scale, below which we can trust the effective field theory, or M = Mp 
(in the case of supergravity). 

In addition, there are also inflaton-induced supergravity corrections to the flat 
direction potential. All these terms provide a general contribution to the flat direction 
potential which are of the form 20 



V{ip) ^ HH4f , Vi^) ^ HMlf , (13) 



^"l , V{^)^HMlf(^ 

where / is some function. The first contribution usually gives rise to a Hubble induced 
correction to the mass of the flat direction, c^fiJ^jcpp, with an unknown coefficient, 
Cjy, which depends on the nature of the Kahler potential. The second contribution is 
the Hubblc-induced A-term. 

Note that ch can have either sign, li ch ~ 1, the flat direction mass is > H. 
It therefore settles at the origin during inflation and remains there. Since \ip\ = at 
all times, the flat direction will have no interesting consequences in this case. The 
positive Hubble induced mass to the flat direction has a common origin to the Hubble 
induced mass correction to the inflaton in supergravity models. This is the well known 
7y-problem |21) . which arises because of the canonical form of the inflaton part of the 
Kahler potential. A large i] generically spoils slow roll inflation. In order to have a 
successful slow roll inflation, one needs t] <^ 1 

The Hubble induced terms can be eliminated completely from inflation and MSSM 
flat direction sectors completely if there is a shift symmetry in the Kahler potcntialsjj, 
or, the Heisenberg symmetry Note that string theory generically gives rise to No- 
scale type Kahler potential [JS] which preserves Heisenberg symmetry at tree level. In 

tt A nice realization of chaotic inflation within supergravity is obtained by implementing a shift 
symmetry I22| . If the inflaton Kahler potential has the form K = (</) -|- </)*)^ /Mp, instead of the 
minimal one K = 4>*ij>/M'^, the scalar potential along the imaginary part of 4> remains flat even for 
Transplanckian fleld values. Therefore it can play the role of the inflaton in a chaotic model. Note 
that a shift symmetry also ensures that the (positive) Hubble induced corrections to the mass of flat 
directions vanishes at the tree-level and the terms in Eq. 1131 disappear. It is also possible to realize 
chaotic inflation for sub Planckian field values in supergravity. For example, see the multi-axions 1231 
driven assisted inflation 1241 . 
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either cases the inflaton and the MSSM sectors are free from Hubble- induced mass and 
A-terms. In the fohowing we wih consider examples where there are no Huhhle-induced 
terms in the potential. 

The relevant part of the scalar potential is then given by 



where A„ ^ and n > 4. Note that the low energy j4-term is a dimensionful 

quantity, i.e., A ~ mg ~ 0(100 GeV — 1 TeV), and depends on a phase. As shown 
in ^SIj all of the MSSM flat directions are lifted by higher-order terms with n < 9. If 
a flat direction is lifted at the superpotential level n, the VEV that it acquires during 
inflation will also depend on the presence of the A-term. 

However the angular direction of the potential being flat, during inflation it 
obtains random fluctuations. There will be equally populated domains of Hubble 
patch size where the phase of the A-term is positive and negative. In either case 
during inflation, the flat direction VEV is given by: 



However there is a distinction between a positive and a negative phase of the A-terms. 
The difference in dynamics arises after the end of inflation. In the case of positive A- 
term the flat direction starts rolling immediately, but in the case of a negative A-term, 
the flat direction remains in a false vacuum for which the VEV is given by Eq. H15() . 
The mass of the flat direction around this false minimum is very small compared to 
the Hubble expansion rate during inflation, i.e. (3ri^ — 9ri + 8)mQ ^ where n > 3. 
During inflation the flat direction obtains quantum fluctuations whose amplitude is 
given by Eq. © . 

The flat direction can exit such a metastable minimum only if thermal corrections 
are taken into account. From now onwards we will assume that the flat direction is 
locked in a false vacuum with a VEV ^pinf during and right after inflation. In the next 
Section we will discuss the various thermal corrections a flat direction obtains. 

5. Thermal corrections to the flat direction 

If part of the SM gauge group remains unbroken by the flat direction VEV, the 
associated gauge fields and gauginos will not receive an induced mass. Together with 
the light particles and sparticles, they reach full thermal equilibrium through gauge 
interactions (for related studies, see refs. |27[ I28| '). The most important processes are 
2 — i- 2 and 2 — + 3 scatterings with gauge boson exchange in the t-channel For 
massless gauge bosons (as happens for the unbroken subgroup) these scatterings are 
extremely efficient and lead to an almost instant thermalization of particles upon their 
production in inflaton decay |2H] tt- 

The presence of a thermal bath with a temperature T affects the flat direction 
dynamics. This happens through the back-reaction of flelds which are coupled to the 
flat direction jllL [T^ . The flat direction VEV naturally induces a mass y(pinf to the 
field which are coupled to it, where ?/ is a gauge or Yukawa coupling. Depending on 
whether yipmf < T or yipinf > T, different situations arise. 

ft However particles carrying quantum numbers associated with the broken gauge subgroup reach full 
equilibrium much later. The reason is that corresponding gauge fields (and gauginos) acquire a mass 
from the flat direction VEV which suppresses the thermalization rate of these particles El 151 . 




(14) 




(15) 
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• y^inf ^ T: Fields which have a mass smaller than temperature are kinematically 
accessible to the thermal bath. They will reach full equilibrium and result in a 
thermal correction Vth to the flat direction potential 

V^th ^ WT'^Wt (16) 

The flat direction then starts oscillating, provided that yT > H ^1 . For iJ < Fd 
the Universe is in a radiation-dominated phase, and hence T oc H^^^. For H > 
the inflaton has not completely decayed and the plasma from partial inflaton 
decay has a temperature T ^ (_ffFdM|)^^^ which implies that T cx H^^'^. 
Therefore in both cases yT > H will remain valid once oscillations start. 

• y^Pinf > T: Fields which have a mass larger than temperature will not be in 
equilibrium with the thermal bath. For this reason they are also decoupled 
from the running of gauge couplings (at finite temperature). This shows up as a 
correction to the free energy of gauge fields, which is equivalent to a logarithmic 
correction to the fiat direction potential 

V^th - ±« THn (|(^|') , (17) 

where a is a gauge fine structure constant. Decoupling of gauge fields (and 
gauginos) results in a positive correction, while decoupling of matter fields (and 
their superpartners) results in a negative sign. The overall sign then depends on 
the relative contribution of decoupled fields f. Obviously only corrections with 
a positive sign can lead to flat direction oscillations around the origin (as we 
require). Oscillations begin when the second derivative of the potential exceeds 
the Hubble rate-squared which, from Eq. IjlTII . reads (^aT^ / (pinf) > H . 

Note that thermal effects of the first type require that fields which have Yukawa 
couplings to the flat direction are in full equilibrium, while those of the second type 
require the gauge fields (and gauginos) be in full equilibrium. If all of the SM group is 
broken by flat direction(s), the reheated plasma will reach full equilibrium much later 
after the inflaton decay [3] • Instead it enters a long phase of quasi-adiabatic evolution 
during which the plasma remains dilute. This implies that fields which are coupled to 
the flat direction, and hence have a large mass, decay quickly and are not produced 
again as the inverse decays are inefficient. As a consequence, thermal effects, which 
crucially depend on the presence of these degrees of freedom, will be weakened. This 
is the reason why we want a subgroup of the SM gauge symmetry to remain unbroken. 

Let us now examine the thermal effects for the curvaton candidates identified 
in Section III. Note that the instantaneous temperature of the thermal bath from 
(partial) inflaton decay is T < for a perturbative inflaton decay. Therefore the 
kinematical condition for inflaton decay via the ^HuHd and ^H^'J' terms, given in 
Eqs. (|10l 111!) respectively, implies that only flelds with a Yukawa coupling A < 10~^ 
to the fiat direction can be in equilibrium with the thermal bath. We only have two 
possibilities: 

• Uidd: S'J7(2)vv remains unbroken in this case. This implies that the correspon- 
ding gauge fields and gauginos, and Hd (plus the Higgsinos) and the LH 
(s)leptons reach full thermal equilibrium. The back-reaction of results in a 

f For example, consider the HuH^ flat direction. This direction induces large masses for the top 
(s)quarks which decouples them from the thermal bath, while not affecting gluons and gluinos. 
Therefore this leads to a positive contribution from the free energy of the gluons. 
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thermal correction AfT-^|(/7|^, see Eq. ®, with Ai ~ 10^^. The free energy of the 

SU (2)w gauge fields result in a thermal correction ~ +awTHii (^{(fl^^- Note that 

the sign is positive since the flat direction induces a mass which is > T (through 
the d) for the LH (s)quarks but not the SU{2)w gauge fields and gauginos. 

• QiL'd: An SU{2) subgroup of SU{3)c is unbroken in this case. Therefore only 
the corresponding gauge fields and gauginos plus some of the (s)quark fields reach 
full thermal equilibrium. Then the back-reaction of Ui and di results in a thermal 
correction (A^ + Aj) T'^l'^l^ according to Eq.©, where A i ~ A2 ~ 10~^. Note that 
logarithmic thermal corrections will not be useful in this case as decoupling of a 
number of gluons (and giuinos) from the running of strong gauge coupling results 
in a negative contribution to the free energy of the unbroken part of SU(3)c- 

Therefore, within MSSM and from the point of view of thermal effects, the Uidd 
flat direction is the most suitable curvaton candidate^ 

6. Uidd as the MSSM curvaton 

The flat directions represented by the udd monomial are partly lifted by n = 4 
superpotential terms, and eventually lifted at the n = 6 level JB]- The former are lifted 
by the udde superpotential term which implies a zero A-teim because udd achieves a 
large VEV which renders a large mass for e and consequently (e) = 0. Therefore the 
scalar potential along these directions will have no metastable minimum. The latter, 
on the other hand, are lifted by the (udd)^ superpotential term which also induces a 
non-zero A-teim (as required). As we shall see, for these directions sufficient thermal 
corrections are induced to lift the flat direction from its false minimum. 

The only channel for two-body inflaton decay is through the $Hu'S' term. The 
existence of two generations of d in the flat direction kinematically blocks two-body 
decay via the other renormalizable term $HuHd (see the discussion in Section O. 
Note that the SU (2)w gauge fields and gauginos are massless as this subgroup remains 
unbroken by the Uidd direction. Therefore degrees of freedom which carry SU{2)yf 
gauge quantum numbers, and their induced mass by the curvaton is < T, fully 
thermalize instantly (compared with the Hubble expansion rate) via 2 — > 2 and 2-^-3 
scatterings with S'C/(2)-w gauge bosons in the t-channel |28l 15). On the other hand, 
the SU{3)c X U{1)y gauge fields and gauginos, plus the right-handed (s)quarks and 
(s)leptons, reach full equilibrium much later as their gauge interactions are suppressed 
by the fiat direction induced mass P). In the following T refers to the temperature of 
SU (2)w degrees of freedom which are in full thermal equilibrium. 

From the discussion in the previous Section it follows that thermal effects result 
in a correction to the effective thermal flat direction potential given by 

^th~ A?TVl'+awr4ln(|¥'l') , (18) 
where Ai ^ 10^'°^ and aw ^ 10^^. 

J As pointed earlier, the L2L3ei fiat direction can obtain a large VEV in the specific case where 
4" = Li. The SU {3)c part of the SM gauge symmetry remains unbroken for this flat direction, and 
hence gluons, giuinos and (s)quarks will reach full equilibrium. We note that neither of L and e 
are coupled to the color degrees of freedom. This implies that there will be no T'^<p^ or logarithmic 
correction to the flat direction potential in this case. This excludes the L2L3ei flat direction from 
being a successful curvaton candidate. 
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Note that the first and foremost condition for fitting the fiat direction from its 
metastable minimum is that Vth > fn-ofinf the other hand, a perturbative inflaton 
decay yields a radiation-dominated Universe whose temperature is T < < m^. The 
reheat temperature Tr ~ \/lj\iMp) is the largest temperature of the Universe in a 
quasi-thermal radiation-dominated phase, where Fd is the inflaton decay rate given in 
Eq. ®. 

If Tr ^ 777,0, the first term on the right-hand side of Eq. (|18|1 dominates over the 
flat direction energy density given by Eq. H14|l for n = 6 right after the inflaton has 
completely decayed. Note that ipinf ^ 10^771-0 from Eq. ©. Note however that very 
soon we wifl have T < Xiipinf since the Hubble expansion rate is gradually decreasing 
and the thermal contribution, i.e. AfT^|(/?|^, becomes ineffective, since the fields with 
coupling Ai to the uidd direction have a mass > T and hence drop out of quasi- 
thermal equilibrium. Eventually the logarithmic thermal correction would take over 
very quickly, even if § Tr « m^. For this reason we focus on the second term of 
Eq. (|18l) in the discussion below. 

First of all we need Vth > iT^afinf' ^'^ ^hat the thermal effects will overcome the 
potential barrier. This leads to the requirement 

awT'^ > mlff^f . (19) 

Second, the thermal mass should trigger flat direction oscillations. If oscillations do 
not begin, the flat direction simply sits at a field value ~ (pinf while Vth is redshifted 
oc in a radiation-dominated Universe. The flat direction will be trapped again in 
the metastable minimum when Vth < "^o'^L/- 

In order for the flat direction oscillations to start we must have cfVth/d\(p\'^ > H^. 
This leads to the condition 

a'4'^>H{T), (20) 

finf 

which always holds in a radiation-dominated phase where H ~ T^/Mp (note that 
Vin/ ^ Mp). This implies that the Uidd direction starts oscillating once the condition 
given in Eq. H19|l is satisfied. This happens, when temperature of the Universe is given 

by 

1/2 

X 10^ GeV. (21) 



{-. 



Vinf 



.1014 GeV^ 

After taking into account of the constraints Tqsc < "m^, which holds for a perturbative 
inflaton decay, and 777^ ~ 10~^ ipinf, which arises from the amplitude of the CMB 
anisotropy, we obtain 

ip^nf > 10" GeV. (22) 

Interestingly enough for directions which are lifted at the n = 6 superpotential level 
we have ipinf ~ 3 x lO^^ GeV, see Eq. (|15|) . Then the infiaton mass is determined from 
Eq. © to be 771^ ~ 3 x 10^ GeV. Also the requirement that Tr > Tosc (a radiation- 
dominated Universe is formed only after the inflaton has completely decayed) restricts 
the inflaton decay rate according to Eqs. © , H21() . A lower bound on the renormalizable 
inflaton coupling is h > 10"* from Eqs. Q), which is perfect for the validity of 
a perturbative infiaton decay Note that the flat direction VEV prevents non- 
perturbative inflaton decay as described in Ref. |14j . 

§ Obviously the T^\ip\^ correction will not arise at all if Tj^ <C m^, since in this case we will have 
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Last but an important reminder to the readers is that the evaporation of the flat 
direction is not an important issue in our case. The interaction term that leads to 
evaporation is of the form Xfip'^x^, with Lp being the flat direction and, x collectively 
denotes the fields in thermal equilibrium. The rate for evaporation is then given by: 



where T^, ^ T and ~ T'^/{v)- Note that the latter being the mass of 

flat direction due to logarithmic thermal corrections, where (</?) is the the amplitude 
of the flat direction oscillations. We then find: 



since Ai ~ 10~^ for uidd flat direction. It takes a large number of oscillations for 
the flat directions to evaporate. Further note that during radiation domination, the 
flat direction VEV scales like: ((p) cx H^/^^ and [ip) cx H during matter domination 
and/or when the flat flat direction oscillates. Therefore, thermal evaporation is not 
an important threat before the curvaton oscillations dominate the Universe ||. 

Hence wc conclude that the udd flat direction is a realistic curvaton candidate in 
the MSSM. At the time of its decay, it dominates the energy density of the universe 
so that its isocurvature perturbations are converted into the observed curvature 
perturbations in CMB. During the last stages of the udd oscillations, the flat direction 
generically fragments into Q-balls [301. The Q-balls eventually decay into LSPs 
through surface evaporation, which can lead to the observed dark matter. Since the 
Q-balls behave like non-relativistic matter, they again start dominating the radiation 
energy density until they completely evaporates. However, this does not modify the 
super-Hubble perturbations. The fragmentation of the flat direction into Q-balls is 
strictly a sub-Hubble process. 

7. Conclusion 

We have considered thermal corrections in the presence of a MSSM flat direction 
condensate, which gives rise to gauge symmetry breaking and a slowing down of 
thermalization rates. We have paid particular attention to the coupling of the inflaton 
to ordinary SM matter and identified among the nearly 300 MSSM flat directions 
the udd direction lifted by n = 6 non-renormalizable terms as a successful MSSM 
curvaton candidate that survives the constraints of energy density dominance and the 
condensate non-dissociation by the ambient plasma. To our knowledge this is the first 
paper where thermal corrections to the MSSM flat direction curvaton are accounted 
for properly. 

The Uidd direction provides masses for the gauge bosons of SU{i)c x C/(1)y 
while leaving S'f/(2)w unbroken. It is thus only the S'[/(2)w degrees of freedom 
that contribute to the thermal correction of the flat direction potential, as given in 
Eq. (|18|l . The dominance of udd over thermal plasma comes about because of two 
important factors. First, the condensate has an initial amplitude which is large, i.e. 
Vm/ ^ 3 X lO^'' GeV. Second, there is an A-term which creates a false vacuum during 
and after inflation which traps the flat direction VEV for a sufficiently long time. 

The CMB fluctuations imparted by udd not only have cosmological implications 
but also astrophysical ones. The flat direction is also a well motivated candidate 

I Along the same lines, it is also true that in the case of a quadratic thermal correction, the 
evaporation rate is again subdominant compared to the Hubble expansion rate. 




(23) 



T^\\{^) «H{T). 



(24) 
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for generating the cold dark matter through Q-ball evaporation. Future coUider- 
based and astrophysical experiments will hopefully pin down the nature of dark 
matter and the physics beyond the SM, but for the time being we are assured that 
minimal supersynimetric Standard Model can provide a cosmologically viable curvaton 
candidate, udd responsible for the CMB fluctuations and reheating, and possibly also 
accounting for the dark matter. 

As in all curvaton scenarios, the spectral index Us for the udd curvaton is very 
close to 1 whereas WMAP 3-year data indicates [2] for the CMB power spectrum 
Hs = O.OSllg Q^g. However, in case of the udd curvaton the spectral index of the 
power spectrum depends also on the yet unspecified inflaton sector with possible tensor 
perturbations. 
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